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Hot carriers in a doped graphene under de electric field is described taking into account the intraband energy 
relaxation due to acoustic phonon scattering and the interband generation-recombination transitions caused 
by thermal radiation. The consideration is performed for the case when the intercarrier scattering effectively 
cstablishes the quasiequilibrium electron-hole distributions, with effective temperature and concentrations of 
carriers. The concentration and energy balance equations are solved taking into account an interplay between 
weak energy relaxation and generation-recombination processes. The nonlinear conductivity is calculated 
for the momentum relaxation caused by the elastic scattering. The current-voltage characteristics, and the 
transition between bipolar and monopolar regimes of conductivity are obtained and analyzed, for different 
temperatures and gate voltages. 
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I. INTRODUCTION 

Study of the heated under de electric field nonequi- 
librium carriers in graphene is stimulated by attempts 
of realization of cfneient ficld-cffect transistor, see 1 - and 
Refs. therein. Essential heating also takes place for the 
cleaning of graphene by strong pulses of applied current^ 
and in recent measurements, 3 conducted in electric fields 
E >1 kV/cm. In addition, excited by ultrafast inter- 
band pumping electron-hole pairs are actively studied 4 . 
In connection with such experiments theoretical studies 
of energy relaxation, both analytical 5 and numcrical 6 , 
were conducted as well as studies of relaxation dynam- 
ics after ultrafast photoexcitation^. Recently, investiga- 
tion of carriers heating under de electric field applied to 
an intrinsic graphene has been given in 8 . Main pecu- 
liarities of the processes of heating in graphene (gapless 
and massless semiconductor, which is described by the 
Weyl- Wallace model^) are related with an interplay be- 
tween quasi-elastic energy relaxation on acoustic phonons 
and generation-recombination processes, that are effec- 
tive nearby the cross point of energy spectrum. In a gated 
graphene, when monopolar regime of transport is real- 
ized, contribution of generation-recombination processes 
is suppressed. Effect of the gate voltage on the processes 
of carriers heating and on current-voltage characteristics 
of a gated graphene (Le., a study of the model of field 
transistor located on graphene sheet of large dimensions) 
have not been considered so far. 

In this paper a heating of carriers in doped graphene 
is analyzed theoretically, for the case when difierence be- 
tween the densities of electrons and holes is controlied 
by a gate voltage. Present treatment is limited by the 
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region of energies below the energy of optical phonon 
{ui ~ 0.2eV), when it is necessary to take into ac- 
count the Joule heating produced by de current, in- 
duced by electric field E, as well as following processes 
of relaxation: a) elastic relaxation of momentum due to 
elastic scattering on structural disorder— , which is the 
most effective mechanism of scattering so that anisotropy 
of carriers distributions is small; b) quasi-elastic relax- 
ation of energy on acoustic phonons; and e) generation- 
recombination processes for interband transitions caused 
by thermal radiation (relaxation of energy and concen- 
trations were analyzed for photoexcitation In—). Here 
we consider the regime when the interparticle Coulomb 
scattering dominates over the processes of relaxation of 
the energies and the concentrations, so that the carri- 
ers are described by quasiequilibrium distributions with 
effective temperature and nonequilibrium concentrations 
of electrons and holes. These parameters are defined from 
the balance equations of the concentrations and the total 
energy, where the areal charge density of graphene sheet, 
given by the difference of the electron and the hole densi- 
ties, is controlied by gate voltage. Besides determinatimi 
of the average energies of carriers and their concentra- 
tions, nonlinear current-voltage characteristics are also 
studied below. 



The analysis performed below is organized as follows. 
The basic equations, describing heating of carriers in 
gated graphene are evaluated in Sec. IL Results of calcu- 
lations, including the parameters of quasiequilibrium dis- 
tributions, nonequilibrium concentrations and current- 
voltage characteristics as a function of applied field, gate 
voltage, and temperature are presented in Sec. III. Dis- 
cussion of the assumptions used and concluding remarks 
are given in Sec. IV. In Appendix the collision integrals 
used in present treatment are specified. 



II. BASIC EQUATIONS 

Our description of the heating of bipolar plasma in 
graphene under de electric field is based on quasiclas- 
sical kinetic equations that take into account the scat- 
tering mechanisms listed above. Here we assume that 
v m S> v cc 3> v qe ,r, where v m , v cc , and v qe (y r ) are the 
relaxation frequencies of momentum, intercarrier scat- 
tering, and energy (concentration) , respectively. Due 
to dominance of the relaxation of momentum, the weak 
anisotropie contributions to the electron (e) and the hole 
(h) distributions are given as^ - 



A/i 



kp 
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(gg ■ p) 

(m) 
PVp 



dfkp \ 
" dp ) 



(1) 



Here k = e (k — h) corresponds the upper (lower) 

sign and the frequency of momentum relaxation Vp n — 
v m p^!(j)l c /h)/h is expressed by means of the character- 
istic velocity v m (that defines efficiency of scattering), 
and of the truncation factor for long-range scattering 10 
ty(x) — x~ 2 exp(— x 2 )I\(x 2 ). Here the correlation length 
of static disorder, l c , and the modified Bessel function 
of the first kind, I\(x), are used. The nonequilibrium 
isotropie distributions fk P are defined by kinetic equa- 
tions 



±eE 



dp 



J qe (f\kp) + Jr(f\ P ) + Jcc(f\kp), (2) 



where an overline stands for the averaging over the angle 
in p-plane. The collision integrals in the right hand side 
of Eq. ^ are given in Appendix. 

Obtained from Eq. ([2]) distributions fk p define concen- 
trations of electrons, n e , and holes, n/j, as follows 



rik = 4 



dp 



{2nh} 



7Jkp: 



(3) 



where the factor 4 takes into account degeneration over 
the spin and the valleys. The current density, I, is given 
as 



4e 



dp 



{2nh) [ 



:V p (A/ e p - Af hp ) , 



(4) 



where v p = vwp/p is the velocity of particle with 
momentum p, written using the characteristic velocity 
vw ~ 10 8 cm/s. The sheet charge is defined by the differ- 
ence of concentrations n e — n^ = An s that is controlied 
by the gate voltage V g applied to the back gate, placed at 
a distance d from the graphene sheet. Considering such 
structure as piane capacitor, we obtain relation between 



An s 



and V g as follows: An s 



eV g /(4ir\e\d), where 



is dielectric Constant of SÌO2 substrate. 

Because of predominance of the intercarrier scattering 
in Eq. ^, the symmetric distributions fk P must satisfy 
the following conditions 

/fep+Ap/fcip'-Ap(l - /fcp)(l - fhip') (5) 

= /fep/feip'(l — fkp+Ap)0- — /fcip'-Ap), 



which nullify the factor in integrand of the intercar- 
rier collision integrai (IA.4p . Hence J cc imposes the 
quasiequlibrium distributions 



Jkp 



exp 



vwP - fJ'k 



1 



(6) 



where T c is the effective temperature of carriers, and /!& 
are the electrochemical potentials of carriers. These three 
parameters are related by requirement for the area charge 
density to be Constant and equal to eAn s , which gives the 
electroneutrality condition in the forni 



eAn, = 4e 



dp 



(2nhy 



"\Jep Jhp) 



(7) 



Further, due to conservation of concentrations of elec- 
trons and holes under the intercarrier scattering pro- 
cesses [see Eqs. (|A.4|> ] we obtain the concentration bai- 
ance equation written through the radiative collision in- 
tegrai (|A.2|) 



dp 



{2nh) 2 



^ [ N 2p/p T ( l ~ fep ~ fhp) ~ fepfhp] = 0. 



(8) 

Here N z = (e z — 1) 1 is the Planck distribution function 
of the equilibrium thermal radiation with temperature T 
and we have introduced the characteristic thermic mo- 
mentum pt — T /v\\r- Point out, to obtain the concen- 
tration balance equation ((SJ we have used that: rik is 
not modified by scattering on acoustic phonons (because 
of smallness of the velocity of sound in comparison with 
Vw) and J r (f\p) is independent of k. 

In a similar way, by summing Eq. ([2]) over p and k with 
the weight VwP, we obtain the energy balance equation: 
P(T C ) — Qj. Here the energy losses term, P(T C ), is 
written as follows: 



P(T C ) = -Av w 



àpp 

(2TTh) 2 



Y, J ie{f\kp)+2J r ( K f\p) 



while the Joule heating term, Qj = I • E, is given as 



(9) 



-Aevw 



dpp 



E 



fdAf ep dAf hp \ 
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Using the integration by parts, one can rewrite Eq. (|10[) 
as Qj = a E 2 , where the nonlinear conductivity a is in- 
troduced by the relation I = crE. Under substitution Eq. 
(1) into (4) we obtain the conductivity 



a = cr 



F p =o 



ir, d 

2] dpF ^ 

o 



(11) 



where F p = f ep + fhp and <jq = (2vw/vd)e 2 /ttH is the lin- 
ear response conductivity for the short-range scattering 
in intrinsic graphene. 

Thus, the parameters of the distributions Eq. ([6]) are 
defined by transcendent equations (jTjl- flTTT) . The nonequi- 
librium concentrations n& and the current density I are 
determined by Eq. ^ and Eq. (QJ, respectively. 



MI. RESULTS 



Below we discuss the solutions of Eqs. (f7l)- (|TTj) as well 
as the concentrations and the current-voltage character- 
istics versus the field E, the gate voltage V g , and the 
temperature T. Calculations are performed at temper- 
ature interval T —11 - 300 K for the graphene on SÌO2 
substrato of width d ~ 3x IO -5 cm and for the correlation 
length l c ~10 nm. 



A. Quasi-equilibrium distribution 

The quasiequilibrium distributions given by Eq. ^ 
are determined from the balance equations, (l7))- ([T0|) . 
where it is convenient to use^the dimensionless momen- 
tum, x = v w p/T c , suchthat f kx = [exp{x— Hh/T )+l]~ l . 
Then the charge neutrality condition (J7J obtains the form 
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The concentration balance equation ([5]) is given as 

e 2x-(^+p, h )/T a _ 1 



(12) 
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and it is not dependent explicitly on the field E or gate 

voltage V g . Multiplying the energy balance equation y \ y and 3 y reS p ec tively. 
P(T C ) = Qj by nh 3 v 3 v /(2v qe T^) we obtain its dimen- 
sionless form P{T C ) = Qe- Where from Eqs. (O and 
(TTU|) it follow the dimensionless energy losses 

P{T c )= 1 -^-± £ / dxx*e*-^fl 



FIG. 1. (Color online) Effective temperature T c /T (a) and 
maximal distribution of electrons / ep =o (b) and holes fhp=o 
(e) versus electric field for T =77 K (1), 150 K (2), and 300 
K (3). Solid, dashed, and dotted curves correspond to V g —0 



2r / dxx 3 f ex f, 



hx 



lo L e 

and the dimensionless Joule heating 

4 



^2x-(^+^ h )/T c _ l 
2xT c /T _ 1 



- 1 ,(14) 



Qe = 



Ve 

2 



/ VWPE \ 
dX I J ex -)- Jfix 



fex=0 + fhx=0 

^(Vcx) 2 



(15) 



Here the characteristic momentum pe is introduced as 
p% = (eEh) 2 /v qe v d and i] c = T c l c /hv w . 

For heavily doped case electrons are degenerated and 
holes are nondegenerated (or vice versa). Then the elec- 
tron quasi- Fermi energy, /i e , is given from Eq. (J12I) as 
fx e ~ yfirn e hvw which is weakly dependent of tempera- 
ture. Further, equation that defines T c (e.g., as function 
of E, for given V g and T) is obtained from Eq. (fT5)) as 



T r ~T 



2* {Helo/hVw) 



(16) 



Point out, the condition of heavily doping (in particular, 
n e ^> Uh) we can rewrite as (/z e /T c ) 2 ^> 1. 

In Figs. la, lb, and lc we plot the effective tem- 
perature T c /T and the maximal distributions, / ep =o and 
fh P =o, for different doping levels at V g —0 V, 1 V, and 3 V 
[notice, in Fig. l(c) the solid curves, pertinent to V g =0 
V, are omitted as they coincide with the solid curves of 
Fig. lb]. It is seen from Fig. la that, for given V g , the 
relative increase of the effective temperature of carriers, 
T c /T, with growing E becomes smaller for larger T. Fig. 
lb shows that the decrease of / ep =o with growing E be- 
comes, at given T, slower for larger V g . From Figs. lc 
and lb it is seen that fhp=o, for given V g , quickly de- 
creases as T grows. For heavily doped case, at V g =10 
V and 20 V, in Figs. 2a and 2b we plot T c /T and f ep=0 
as functions of E. For these gate voltages, the hole con- 
centrations are small, so that fhp=o < 0.1 and we do not 
plot the dependencies fh P here. The dependencies T c /T 
and f ep =o on T, V g , and E now are similar to those in 
the low-doping region. 



B. Carrier concentrations 

Using solutions of the balance Eqs. (JT2J - (JT5J) in Eq. 
(|3"|) we plot in Fig. 3a dimensionless concentrations of 
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FIG. 2. (Color online) Eflective temperature T c /T (a) and 
maximal distribution of electrons / ep= o (b) versus electric 
field for T =77 K (1), 150 K (2), and 300 K (3). Solid and 
dashed curves correspond to V g =10 V and 20 V, respectively. 
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FIG. 3. (Color online) Electron (a) and hole (b) concentra- 
tions (normalized to their values at E = 0) versus electric 
field at T =77 K (1), 150 K (2), and 300 K (3). The solid, 
the dash, the dot, the dash-dot, and the dash-dot-dot curves 
correspond to V a =0 V, 1 V, 3 V, 10 V, and 20 V, respectively. 



electrons, n e /n ei E=0i and in Fig. 3b of holes, rih/nh.E=o- 
Here n e ^E=o and Uh,E=o are equilibrium concentration of 
electrons and holes for E = and these concentrations 
are strongly dependent on V g and T*£ For small concen- 
trations (solid and dashed curves in Fig. 3a) of electrons 
(the main carriers) their concentration is increased on 
tens of percents for growing E. However, for V g > 10V 
the concentration of main carriers becomes very weakly 



growing function of E even at T = 300K. At the sanie 
time, the density of the minor carriers (holes) according 
to Fig. 3b can be enlarged by tens of times, especially at 
low temperatures. 
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FIG. 4. (Color online) Nonlinear conductivity, normalized to 
its value at E — 0, versus electric field for: T =300 K (a), 
150 K (b) and 77 K (e). The solid, the dash, the dot, the 
dash-dot, and the dash-dot-dot curves correspond to V g = 
V, 1 V, 3 V, 10 V, and 20 V, respectively. 



C. Current-voltage characteristics 

In Fig. 4, we plot the normalized conductivity a/aE=o 
versus field E for different levels of doping, dependent on 
gate voltage varied between V g — and V g =20 V and 
on different temperatures. Point out &e=o for each curve 
in Figs. 4 a-c has specific value and <7e=o/&o increases 
essentially when V g grows J^ The field dependency of the 
conductivity is determined by a competition between the 
first and the second terms in Eq. (fTTj) . Indeed, as F p =o 
is a decreasing function of E the integrai (second) contri- 
bution, for [i e l c /hv\v <] 1 or/and T c l c /hvw ^ 1, is either 
much slower decreasing function of E or even an increas- 
ing function of E. In particular, for T c l c /hvw > 1 and 
fi e /T c ^> 1, the second contribution is the integrai of the 
product of F p (that is « -F p =o for /J, e /vw > P > and 
drops rapidly for p > fi e /vw) with the other factor that 
is a large and growing function of p, oc p 2 , in actual re- 
gion of p, T c /vw ;$ P < He/vw- So for low temperatures 



the contribution F p=0 (see Fig. 4 b and e) is essential 
and leads to the appearance of a minimum in the field 
dependence of a for small concentrations, for V g < 3V. 
With the growth of electric field or the concentration 
(for V g >3 V) the main contribution to the conductivity 
comes from the integrai term in Eq. (|11|) and now a is 
monotinically growing function of E, cf. Fig. 4 a, b, 
and e. Indeed, pertinent current-voltage characteristics 
show transition from a sublinear dependence to a super- 
linear one, as the concentration or/and the temperature 
grow. For V g >3 V (and at T =300 K for any V g , or con- 
centration) a superlinear current-voltage characteristic is 
rcalizcd. 



are imposed also because of the lack of data on graphene. 
To conclude, essential heating of carriers is shown in 
present study. It will define different parameters of pos- 
sible graphene-based devices. In addition, investigation 
of hot carriers (even within the considered limited re- 
gion of parameters) gives important information about 
mechanisms of relaxation and recombination in graphene. 
Therefore, further experimental and theoretical (includ- 
ing numerical modeling) study of a heating in graphene 
is very opportunely now. 



Appendix: Collision integrals 



IV. CONCLUSIONS 

In the present work, we investigate the effect of doping 
on carrier heating in a gated graphene based on the bai- 
ance equation approach. The effective temperature and 
concentrations of carriers are studied as functions of the 
thermostat temperature, the gate voltage, and the ap- 
plied de field. Pertinent current-voltage characteristics 
show the transition from the regime of heating for intrin- 
sic bipolar conductivity to the regime of heating for the 
conductivity of strongly degenerated carriers. In the lat- 
ter case the heating is well described by simple analytical 
formulas, see (16). 

Next, we list and discuss the assumptions used in 
the calculations performed. Here we have examined 
the heating of carriers in the low energy region only, 
when the optical phonon emission is not essential and 
the radiative-induced direct interband transitions are as- 
sumed to be the main generation-recombination mech- 
anism. Point out that for much larger electric fields 
(for E 3>0.1 kV/cm, see^) present approach will not 
be valid as the optical phonon contribution is not taken 
into account, however, the influence of optical phonons 
is smaller for a higher level of the doping (gate voltage) . 
We have also restricted ourselves by the study of limiting 
case when the intercarrier Coulomb scattering is domi- 
nating. As shown in 8 for the intrinsic graphene case, 
this is an adequate approach. Possible contribution of 
other generation-recombination mechanisms (note, that 
the Auger-processes are forbidden due to the symmetry 
of electron-hole statesi) require an additional investiga- 
tion. Further, we have assumed that the heat removal 
is sufficicntly effective (this point becomes insignificant 
in the case of short electric pulses when the thermostat 
is not overheated). At last, to describe the momentum 
relaxation we take into account only the sfaticai disor- 
der scattering using the phenomenological model of Ref. 
10 (which is in good agreement with the experimental 
data; the microscopie mechanisms of scattering are stili 
unclear— ). The listed assumptions should not essentially 
modify the peculiarities of the heating mechanisms, in 
comparison with a study that will treat relaxation pro- 
cesses in more details; in addition, here many limitations 



Below we present the collision integrals J qe , J r and J cc 
(deduced i n 10 i n andì&, respectively) that have been used 
to obtain the balance equations. For quasi-elastic scat- 
tering on acoustic phonons the energy relaxation is de- 
scribed by the Fokker-Planck nonlinear differential form 



Jqe (f\kp) 



,(?e) 



p 2 dp 



dfkp 
dp 



/fcp(l _ fkp) 



Pt 



(Al) 
Here k = e,h and it is introduced the frequency of quasi- 



elastic relaxation u. 



V 



'q e p/h that contains the char- 
acteristic velocity v qe oc T, where for typical parameters 
of graphene we have v qe s=a 51.4 cm/s at T = 300K. The 
thermal-radiation-induced interband transitions are be- 
ing depicted by the collision integrals 

Jr (f\p) = 4 T) [ iV 2 P / PT (l ~ fep ~ fhp) - fepfhp] , (A.2) 

that are equal for electrons and holes. In Eq. (|A.2|) 

(r) 

it is used the frequency of radiative relaxation v v = 
v r p/h, where the characteristic velocity v r ~ 41.6cm/s 
for graphene imbedded between SÌO2 substrate and cover 
layer. 

As Auger-processes are forbidden because of the 
electron-hole symmetry of the bands in graphene^, the 
intercarrier collision integrai receives the form 

00 00 

Jcc(f\k P )= J2 J dp' J dA P W kkl (pp'\Ap) 

ki=k,k' q _ 0o 

X [fkp+Apfkip'-ApO- - /fcp)(l _ /fcip')(A-3) 
— /fcp/feip'(l — /fcp+Ap)(l — fkip'-Ap)] , 

where it describes the transitions due to carrier-carrier 
scattering from the initial states (fep, fc'p') to the states 
(fcp + Ap, fc'p' — Ap). Here uvHApI determines transfer- 
enee of the energy for scattering and k = k' corresponds 
to electron-electron or holc-hole transitions whereas the 
channel k ^ fc' describes scattering of electrons on holes. 
Also in Eq. (|A.4[) it is carried out average over the angle 
and as a result the probability of transition Wkk' (pp'ì^p) 
is dependent from |p|, |p'| and from the transference of 
the energy vyyAp. Deduction of the concentration ([8|) 
and the energy ((5|) balance equations, needed, along with 



Eq. ©, to obtain the parameters of quasiequilibrium 
distributions Eq. ©, is based on the property of con- 
servation as of the concentration so of the energy density 
for interparticle scattering as 



4 

L 2 



p 



4vw 
1? 



J2pJ cc (f\kp)=Q. (A.4) 



These conditions straightforwardly follow from Eq. (|A.4j) 
if take into account the symmetry of the probability of 
transitions under the interchange of scattering channels 
k, p and k', p'. 
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